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1.0  Introduction 


1 . 1  Purpose  and  Scope 


The  objective  of  this  effort  is  to  experimentally  measure  the 
current  and  charge  distributions  of  various  analytical  models,  such 
as  cylinders,  crossed  cylinders  and  a  scale  model  747  aircraft,  in 
free  space  and  in  proximity  to  a  perfectly  conducting  ground  plane. 
Such  experimental  data  will  aid  in  better  understandinn  the  mutual 
interaction  of  the  Horizontally  Polarized  Dipole  and  Vertically 
Polarized  Dipole  (HPP/VPD)  simulators  and  test  objacts  as  well  as 
providing  a  data  base  for  comoarison  with  calculations. 

The  report  describes  results  obtained  from  a  series  of  time-domain 
scale  model  electromagnetic  measurements.  The  tests  were  conducted 
at  the  Lawrence  Livermore  National  Laboratory  (LLNL)  Transient 
Electromagnetic  Range  and  involve  the  use  of  subnanosecond 
electromagnetic  pulses  to  illuminate  the  bodies  under  test.  Small  B 
and  5  sensors  were  used  to  measure  the  local  fields  on  the  bodies  to 
determine  the  induced  currents  and  charges. 

1 .2  Background 

Major  concerns  in  the  evaluation  of  a  full-scale  simulator  such  as 
the  Air  Force  Weapons  Laboratory  (AFWL)  VPD  and  HPD  include  more 
accurate  knowledge  of  simulator-test  object  interaction  ahd  the 
relationship  between  the  response  of  the  object  in  the  simulator  and 
its  response  in  free  space.  The  full-scale  simulator  data  are 
extrapolated  to  object  response  in  the  actual  operational 
environment  by  use  of  extrapolation  functions.  These  functions  are 
often  generated  b^  means  of  numerical  models  which  predict  response 
of  objects  in  free  space  (if  such  is  the  operational  environment) 
and  the  response  when  the  object  is  in  close  proximity  to  a  ground 
screen  (which  is  how  the  ground  of  the  VPD  constructed).  An  onqoinq 
need  therefore  exists  for  experimental  data  to  confirm  the  validity 
of  numerical  models. 


One  method  for  qeneratinq  such  data  is  with  transient  scale-model 
measurements.  It  is  only  within  the  past  few  years  that 
subnanosecond  Dulse  qeneration  and  sampling  technology  has  advanced 
to  the  point  where  useful  scale  model  time-domain  electromagnetic 
pulse  (EMP)  measurements  may  be  performed  in  the  laboratory.  The 
Lawrence  Livermore  National  Laboratory  has  such  a  facility  which  was 
used  to  obtain  the  data  presented  in  this  report. 

1 .3  Organization  of  Report 

In  Section  2,  data  collection  techniques  are  described.  Included  is 
a  description  of  the  operation  of  the  transient  range,  description 
of  sensors  and  data  collection  procedures.  In  Section  3,  the  test 
objects  in  question  are  described.  These  include  a  cylinder, 
crossed  cylinder  and  the  747  aircraft  model.  In  Section  4,  critical 
issues  in  the  collection  of  the  raw  data  and  the  encountered 
difficulties  are  described.  Data  nrocessing  techniques  follow  in 
Section  S.  The  data  for  the  three  test  objects  (both  raw  and 
processed)  are  summarized  in  tabular  form  in  Section  6.  The  bulk  of 
the  data  is  relegated  to  apoendices.  Comparisons  with  test  data 
from  the  University  of  Michigan  and  with  numerical  prediction  follow 
in  Section  7.  Results  for  pole  extraction  are  summarized  in 
Section  8. 

2.0  Data  Collection  Techniques 

2.1  LLL  Transient  Range 

The  LLNL  transient  ranqe  facility  [Ref.  1]  consists  of  a  monocone 
antenna  that  radiates  electromagnetic  pulses  over  a  horizontal 
aluminum  ground  plane.  A  photograph  of  the  range  is  shown  in  figure  1. 


1.  Oeadrick,  F.  Miller,  E.  K.,  and  Hudson,  H .  G.,  The  LLL  Transient- 
Electromagnetics-Measurement  Faci 1 ity,  Lawrence  Livermore  National 
Laboratory,  Livermore,  CA,  Rept.  UCRL-51933  (197S). 
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The  ground  plane  is  8.5  m  x  8.5  m  large.  Radiated  electric  and 
magnetic  field  measured  on  the  ground  at  ? . 4 3  m  from  the  base  of  the 
monocone  antenna  are  shown  in  Figures  2  and  3  (this  corresponds  to 
the  location  of  test  objects  where  the  experimental  data  are 
obtained).  The  rise  time  of  the  pulse  is  topically  less  than  100  ps 
and  the  pulse  width  is  about  300  ps.  Reflections  from  the  edges 
have  been  found  to  be  small  and  the  ranqe  has  been  used  with  data 
records  up  to  200  ns. 

Additional  information  pertinent  to  the  data  reported  here  is  given 
in  Appendix  A,  which  also  includes  a  description  of  range 
improvements  that  were  undertaken  in  conjunction  with  this  effort. 
Validation  and  calibration  issues  concerning  the  range  have  been 
dealt  with  elsewhere  [Ref.  2j. 

2.2  Field  Sensors:  Roth  tanqential  maanetic  fields  and  normal  electric 
fields  on  the  test  objects  were  measured.  These  were  performed  with 
the  AC D- 1  D  sensor  (Ref.  3)  and  the  MGL-8  B  sensors  (Ref.  A).  3oth 

of  these  precision  sensors  have  been  precalibrated;  the  calibration 

_  a  ■?  .  -  5  9 

factor  is  Ap^  -  10  m"  for  the  D  sensor  and  Apq  -  10  m“  for 

the  8  sensor. 

Figures  2  and  3,  which  show  the  incident  E  and  H  fields  at  the  test 
location  on  the  ground  plane,  were  determined  from  measurements  with 


2.  Bevensee,  R.  M.,  Deadrick,  F.  J. ,  Miller,  E.  X.,  and  Ok  ad a,  3.  T. ,  VaJ ;dat >Q' 
and  Calibration  of  the  LLL  Transient  Electromagnetic  Measurement  ~a:i'uv, 
Lawrence  Livermore  National  Laboratory,  Livermore,  CA,  Reot.  HCRL-5222-i  (19771. 

3.  Baum,  C.  E.,  Breen,  E.  L.,  Giles,  .1.  C.,  O’Neill,  J.,  and  Sower,  G.  0., 

"Sensors  for  Electromagnetic  Pulse  Measurements  Both  Inside  and  Away  from 
Nuclear  Source  Regions,"  IEEE  Trans.  Antennas  Propaq.  AR-26,  No.  1,  January, 
1978. 

4.  Olsen,  S.  L.,  "Sensor  MGL-98  sensor  OW,"  AL-1186,  September  1975,  Albuquerque 
Division,  EG&G  Inc.,  9733  Coors,  R.,  N.  R . ,  Albuquerque,  NM. 
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Figure  ].  LLNL  Transient  Antenna  Range  showing  the  Ground  Plane  and  the 
Monocone  Source  Antenna. 
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these  sensors.  The  ibsolut'1  eloctri'-  r  hi  is  found  by  nutner ical  1  y 
integrating  the  v'w  wutput  of  the  0  sensor: 


£(t) 


e  A  z. 
o  eq  I. 


(tl  dt  , 


(1  ! 


where  E(t)  is  the  electric  field  in  V/m  os  a  function  of  time, 

is  the  permittivity  of  free  spa.e,  z,  is  the  impedance  of  the  load 

connected  to  the  sensor  output,  (nominally  50  here),  and  A  is 

•  * 
the  equivalent  area  of  the  D  sensor.  Similarly,  the  magnetic  field 

can  be  found  from  the  output  of  the  5  sensor  by  the  expression 


H  ( t )  = 


H-  A 
o  eq 


i: 


V  „  (t)  dt 

out*  1 


t  ?  ■ 


where  H(t)  is  the  magnetic  field  in  amperes/meter ,  p  is  the 

permeability  of  free  space,  A  is  the  equivalent  area  of  the  5 

^  • 

sensor  and  VQut.(t)  the  output  of  the  5  sensor  as  a  function 
of  time. 


Particularly  qood  agreement  in  both  maqnitude  and  temooral  shape  is 
obtained  when  we  inteqrate  and  calibrate  the  sensor  outputs.  ('Jote 
that  in  free  space  F  -  q  H ,  where  q  is  the  imoedance  of  free  soacel 
Typicallv,  the  peak  electric  field  at  the  test  ooint  is  on  the  order 
of  ABO  V/m  and  the  peak  magnetic  field  >  1.3  A/m,  when  the  conical 
antenna  is  driven  with  an  impulse  generator  with  a  peak  voltage  of 
1.5  kV. 

•  • 

A  Dhotoqraoh  of  the  D  and  B  sensors  used  is  shown  in  Figure  4. 


3.0  Description  of  Test  Objects  and  Experimental  Setup 


In  these  experiments,  two  types  of  measured  electromagnetic  quantities 
were  obtained:  the  normal  electric  field  and  the  surface  maqnetic  fields 
at  soecific  points  on  the  scale  models.  These  measurements  were  obtained 


Figure  2. 


Figure  3 
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Measured  Electric  Field  at  Test  Point  2.43  m  from  Conical  Antenna 
on  the  Image  Plane. 
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for  objects  in  both  a  "free-space"  configuration  and  in  the  presence  of  a 
perfectly  conducting  qround. 

The  free-space  measurements  are  actually  a  simulated  free  space  where  a 
plane  of  symmetry  is  used  as  shown  in  Fiqure  5(a).  This  configuration 
allows  only  ant i symmetr ic  excitation  on  the  model,  and  thus  is  a  soecial 
case  of  free-space  operation. 

For  the  measurements  of  models  near  a  perfectly  conducting  ground  plane, 
a  second  vertically  oriented  plane  was  used  as  shown  in  Fiqure  5(b).  The 
vertical  qround  plane  was  constructed  so  that  it  could  easily  be  moved  in 
and  out.  For  these  experiments  the  model  remained  fixed  relative  to  the 
source  monocone  antenna;  the  ground  plane  was  moved  as  required.  A 
sketch  of  the  experimental  layout  with  dimensions  is  shown  in  Fiqure  6. 
The  vertical  qround  plane  was  made  of  a  thick  aluminum  plate  with 
dimensions  of  1.2  X  4.  m. 

3.1  Cylinder 

Fiqure  7  is  a  photograph  of  the  cylinder  model  used  in  all  these 
tests.  The  lenqth-to-radius  ratio  for  the  cvlinder  model  is 
L/a  =  20,  where  L  is  the  total  length  of  the  cylinder,  and  a  is  the 
radius.  While  L  is  defined  as  the  total  lenqth  of  the  cylinder 
being  simulated,  only  one  half  appears  above  the  symmetry  olane  used 
in  these  experiments.  The  cylinder  model  was  constructed  of  a 
section  of  brass  pipe  with  a  wall  thickness  of  0.535  cm. 

A  flat  end  cap  was  machined  to  fit  the  end  of  the  cylinder,  and  the 
metal  to  metal  joint  was  sealed  with  a  conductive  silver  compound  to 
insure  a  qood  electrical  junction.  Two  small  flat  surfaces  were 
milled  on  the  surface  of  the  cylinder  to  accommodate  the  R  and  0 
sensors . 

To  obtain  a  good  electrical  contact  between  the  cylinder  and  the 
conductive  symmetry  plane,  a  s i 1 ver- loaded  conductive  paste  was  also 
applied  to  the  bottom  of  the  cylinder  over  the  area  that  comes  in 

1 1 


•  • 

Figure  4.  AFWL  MGL-8  B  and  ACD-4  L)  sensors  used  in  the  time-domain 
measurements . 
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Figure  6.  Experimental  configuration  for  measurements  of  bodies  nea>- 
perfectly  conducting  ground  plane.  (The  distance  h  was  set 
at  a,  5a,  10a,  and  20a  (10  cm,  50  cm,  1  m,  and  2  m].) 


contact  with  the  aluminum  olane.  The  miniature  field  sensors  were 
held  in  position  with  a  cooper  foil  tape  with  a  conductive 
adhesive.  The  dimensions  of  the  cylinder  and  sensor  locations  are 
shown  in  a  s<etch  in  Fiqure  H.  A  top  view  of  the  cylinder  in  the 
ranqe  with  dimensions  is  shown  in  a  sketch  in  Fiqure  1. 

The  rationale  for  usinq  the  half-lenqth  shown  in  Fiqure  3,  which  has 
to  do  with  several  considerat ions  includinq  phase  anqle  errors,  is 
qiven  in  Aooendix  -i . 

J . ?  Crossed  Cy 1 inder 

The  second  model  to  on  tested  is  the  crossed  cylinder  shown  in  the 
photograph  if  Figure  10.  The  dimensions  of  the  model  are  shown  on 
fne  sketch  in  Fiqure  11.  These  dimensions  correspond  to  tnose  of 
the  cylinder  model  used  in  the  previous  section.  Semicircular  end 
caps  were  orazed  to  the  ends  of  the  bottom  section  of  the  cross,  and 
a  flat,  removable  end  cap  was  used  at  the  top  of  the  cross  so  as  to 
allow  connections  to  he  made  to  the  *  sensor  located  at  the  too  of 
the  cross.  The  d  sensor  was  oriented  so  as  to  measure  the  surface 
current  flowing  down  the  vertical  section  of  the  cross;  it  was 
mounted  near  the  junction  of  the  cross  members.  To  ensure  a  nood 
contact  between  the  cross  and  symmetry  olane,  a  thin  sheet  of  brass 
plate  was  epoxied  to  the  base  of  the  cross  with  conductive  epoxy 
(tne  joint  was  also  painted  over  with  silver  paint),  and  the  brass 
plate  was  then  taped  to  the  symmetry  olane  with  conductive  taoe. 

3 . i  /A/  Scale-Model  Aircraft 

A  tnird  series  of  measurements  was  oerformed  on  a  I :  Hid  scale  model 
of  a  doeing  74/  aircraft.  The  model  used  is  a  commercially 
available  elastic  xit  that  is  conveniently  sectioned  into  riqht  anu 
left  nalves,  oP-.iii  lono,  as  shown  in  the  pnotoqranh  of  p i pure  1 
dne-half  of  the  plane  was  assembled  in  a  wtieels-np  configuration  and 
epoxied  to  a  tnm  sheet  of  brass  plate  as  shown  in  the  pnotogriph. 
cour  sensor  mount  inq  holes,  i.'ireo  on  the  wing  top  and  a  fnurtn  no 


Figure  9.  Top  view  of  configuration  used  for  cylinder 
free-space  measurements. 


Removable  end  cap 


Figure  11  ( 


Figure  1 


a).  Dimensions  of  crossed  cylinder  model  used  in  transient 
measurement  tests . 


Perfectly  conducting  ground  plane 


1  (b).  Top  view  of  crossed  cylinder  showing  direction  of  incident 
electric  field. 


the  fuselage  between  the  front  and  rear  winqs,  were  formed  with 

•  • 

expoxy  and  milled  to  form  a  flat  surface  to  accommodate  the  0  and  3 
sensors.  The  entire  model  was  then  qiven  several  coats  of  a 
conductive  silver  paint  to  form  a  thin  metallic  surface.  The 
pertinent  dimensions  of  the  model  are  shown  in  Figure  13. 

4.0  Critical  Issues  in  Haw  Data  Collection 

Much  of  the  data  collected  and  processed  as  part  of  this  effort  is 
characterized  by  a  rather  high  Q.  For  data  processing  purposes,  this  nas 
necessitated  the  taking  of  data  records  over  time  durations  far  in  excess 
of  those  previously  used.  For  examole,  prior  to  this  effort,  the  range 
had  been  used  for  taking  data  records  of  up  to  typically  20  ns.  Some  of 
the  data  collected  and  processed  here  was  recorded  over  a  200  ns 
duration.  This  gave  rise  to  some  problems  that  had  to  be  resolved. 

The  first  area  of  concern  was  edge  reflections  from  the  extremities  of 
the  range.  From  purely  dimensional  considerations,  the  clear  time  was 
established  to  be  20  ns.  However,  data  for  the  simulator  free-field 
taken  over  longer  records  indicated  reflections  far  lower  than  1 0%  of  the 
simulator  output  signal.  The  uncertainty  this  introduced  to  the  data  was 
not  considered  to  be  deleterious. 

The  other  oroblem  encountered  had  to  do  with  drift  in  the  signal 
baseline.  The  data  taken  with  the  combination  of  the  samplinq 
oscilloscope  as  shown  in  Figure  A1  exhibited  excessive  low  frequency 
drift  for  time  durations  exceeding  50  ns.  The  drift  can  introduce 
comol ications  when  processing  the  data.  The  use  of  a  Tektronix  7912 
transient  digitizer  in  combination  with  an  LSI -11  microcomputer  neloed 
eliminate  the  problem.  In  this  report,  data  for  the  objects  in 
free-space  (where  20  ns  records  are  sufficient)  wpre  taken  with  the  older 
system  as  shown  in  Figure  Al.  Longer  data  records  have  all  been  taken 
with  the  Tektronix  7912  digitizer. 

Both  the  older  and  the  new  data  acquisition  system  have  the  potential  of 
sampling  the  data  at  512  points.  Depending  on  the  highest  frequency 
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content  of  the  signal,  this  can  limit  the  total  tune  duration  over  whisn 
tne  data  are  sampled.  For  example,  the  outnut  free-field  of  the 
simulator  could  only  be  sampled  adequately  up  to  a  total  time  duration  j 
50  ns. 

Above  this  a  noticeable  prop  in  peak  amplitudes  results  due  to 
insufficient  time  resolution.  The  object  responses  on  the  other  nano 
were  not  characterized  by  rise  times  as  fast  as  those  of  the  incident 
field  and  no  samplinq  problems  were  encountered  even  for  200  ns  time 
durat i ons. 

The  time  durations  used  for  collection  the  data  varied  doping inc  an 
configuration.  These  are  listed  in  Table  1.  These  time  durations  we-« 
found  to  be  adequate  to  perform  the  required  signal  orocessioq. 

5.0  llata  Processing  Techniques 

Two  sorts  of  siqnal  processing  techniques  were  used.  In  the  first  place, 
the  Fourier  transforms  were  computed  for  all  responses,  and  from  this 
information  the  transfer  function  of  the  objects  (i.e.,  ratio  of  object 
response  to  incident  field)  was  evaluated.  The  time  records  used  to 
obtain  the  frequency  spectra  are  listed  in  Table  1.  For  tne  objects  in 
free  space,  the  data  were  colleited  with  the  older  system  shown  in 
Figure  A1  processed  on  a  CDC-7600  using  - i  Ion’s  Method  (the  Fourier 
integral  is  performed  using  straight  line  interpol ation  to  the  sampled 
time  amplitudes).  The  data  with  longer  time  records  were  recorded  with 
tne  Tektronix  7912  transient  digitizer  in  combination  with  an  LSI-11 
microcomputer.  In  processing  the  data,  a  fast  Fourier  transform 
algorithm  in  the  LSI-11  system  was  used. 

Sesides  tne  freguency  spectra,  the  poles  were  extracted  for  some  of  the 
responses.  The  tecnniques  used  to  extract  tne  poles  are  described  in 
Appendix  C.  Tnree  techniques  were  used;  namely,  irony's  technique,  the 
extended  Kalman  filter  and  the  maximum  likelihood  identifier.  Th« 
trade-offs  between  these  various  technigues  are  discussed  in  great  ootail 
in  the  Apoendix. 


TABLE  1  OBJECT  CONFIGURATION  AND  CORRESPONDING 
DATA  RECORD  TIME  DURATION 


Conf  iquratiuns 


T i me  Duration  ( ns ) 


Simulator  free-field  20 
Object  in  free-soace  20 
Object  2  m  from  qround  plane  50 
Object  1  m  from  ground  olane  50 
Object  50  cm  from  ground  olane  100 
Object  10  cm  from  ground  olane  200 


In  tne  cases  where  time  domain  record  lenqtns  are  Ipss  trian  or  equal  to 
50  ns,  the  processing  is  relatively  simple.  The  time  domain  responses 
are  multiplied  by  a  Hamming  window  function  to  reduce  spectral  leakaqe. 
Spectral  ieakaqe  is  caused  by  tne  calculation  of  the  Fourier  transform 
for  a  finite  data  sequence.  The  Hamming  window  equation  is: 

w  (t)  =  .54  +  .46  cos  (2irt/T  ) 

where  t  is  the  length  of  tne  signal  to  be  windowed.  These  then  are 
transformed  into  the  frequency  domain  and  divided  Oy  the  incident  nulso 
frequency  response.  Tne  100  ns  records  were  processed  in  the  same  way 
except  the  50  ns  incident  pulse  was  used  for  the  division.  The  proper 
frequency  response  for  tne  50  ns  record  can  Do  obtained  by  either  of  two 
methods.  One  method  is  to  straight  line  interpolate  new  frequency  points 
between  the  calculated  points.  The  other  method  is  to  append  zeroes  to 
the  50  ns  record  until  it  is  the  desired  lenqth  and  then  transform  into 
the  frequency  domain.  Eitner  method  can  be  used  to  obtain  the  proper 
frequency  spaced  ooints. 

The  incident  pulse  record  could  be  obtained  only  for  a  50  ns  duration. 
Since  the  transient  dinitizer  collects  a  fixed  number  of  points,  as  the 
record  lenqth  increases  the  time  resolution  bptween  points  decreases. 

This  is  a  problem  for  tne  incident  pulse  because  it  has  very  hiqn 
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frequency  components  arid  a  low  siqnal  amplitude.  Tne  model  responses, 
however,  are  higher  i i .  siqnal  amolitude  <1 nd  do  not  contain  the  high 
frequency  components  so  a  large  time  resolution  is  acceptable. 

6.0  Results 

For  each  of  the  test  objects,  the  data  at  appropriate  positions  were 
collected  with  the  3  and  D  sensors.  The  ambient  E-  and  H-fields  at  tne 
location  of  tne  test  objects  were  also  obtained.  Tne  transfer  functions 
in  the  frequency  domain  were  obtained  by  dividing  tne  test  objects' 
responses  by  the  incident  field  in  tne  frequency  domain. 

The  results  of  these  tests  are  tabulated  in  this  section  and  the  bulk  of 
the  data  is  shown  in  appropriate  appendices.  In  all  the  taDles,  the 
entries  denote  the  following: 

•  Separation  distance  h:  Distance  from  cylinder  and  cross-cvl inder  to 

ground  plane.  See  Figures  5,  6  and  11. 

•  Angle  0  :  Incidence  Angle  in  Degrees.  See  cigure  9. 

’  Peak  :  Peak  amplitude  of  the  transfer  function  at  the  first 

harmonic. 

"  Fig.  #  ;  Indicates  figure  number  for  case  under  consideration 

in  the  appropriate  appendix. 

6.1  Cyl inder: 

The  experimental  predictions  for  the  cylinder  are  tabulated  in 
Table  2. 

6.2  Cross  Cylinder: 


The  experimental  predictions  for  the  cross-cylinder  are  tabulated  in 
Table  3. 


TABLE  2.  SUMMARY  OF  RESULTS  FOR  CYLINDER 


Separation  Angle  .  .  .  transfer  Function 


Distance 

h. 

e 

Oegrees 

D-Sensor 
Fig.  # 

3-Sensor 
Fig.  # 

E/E 1 nc 

Fig.  #  Peak 

H/H 
Fig.  # 

inc 

Peak 

F 

0 

Dl 

1)1 

D6 

9 

Oil 

22 

R 

45 

02 

02 

D7 

9 

D12 

2  ? 

E 

E 

90 

D3 

03 

08 

10 

Dl  3 

20 

S 

P 

135 

D4 

04 

D9 

9 

Dl  4 

20 

A 

C 

130 

D5 

ns 

DIO 

9 

015 

20 

E 

10  •  cm 

0 

D16 

024 

032 

35 

040 

43 

180 

020 

028 

036 

23 

044 

190 

50  •  cm 

0 

017 

D25 

D33 

50 

D41 

56 

180 

D21 

029 

037 

43 

045 

50 

1  •  m 

0 

018 

026 

034 

25 

042 

33 

180 

D22 

030 

038 

24 

046 

27 

2  •  m 

0 

019 

D27 

035 

38 

043 

48 

180 

D23 

031 

039 

37 

047 

45 

TABLE  3.  SUMMARY  OF  TIME  DOMAIN  RESULTS  FOR  CROSS-CYLINDER 


Separation  Angle  .  .  Transfer  Function 

Distance  9  D-Sensor  3-Sensor  t/Einc  H/Hinc 


h. 

Degrees 

Fig.  # 

Fig.  # 

Fig.  # 

Peak 

Fig.  # 

Peak 

Free 

0 

El 

E3 

E5 

23 

E7 

35 

Space 

180 

E2 

E4 

E6 

24 

E8 

28 

17**  cm 

0 

E9 

E 1 7 

E25 

70 

E33 

93 

180 

E13 

E21 

E29 

58 

E37 

130 

50  •  cm 

0 

E10 

E 18 

E26 

46 

E34 

40 

180 

E14 

E22 

E30 

45 

E38 

36 

1  •  m 

0 

Ell 

E19 

E27 

25 

E35 

24 

180 

E15 

E23 

E31 

22 

E39 

16 

2  •  m 

0 

E12 

E20 

E28 

37 

E36 

32 

180 

E16 

E24 

E32 

37 

E40 

29 

*The  cross-cy 1 inder  could  be  brought  only  as  close  as  1 7  cm  to  the  ground  due 
to  physical  obstructions. 


6.3  747  Scale-Model  Aircraft 


The  experimental  predictions  for  the  747  scale-model  aircraft  are 
tabulated  in  Table  4.  The  positions  on  the  model  are  identified  in 
Figure  13.  The  entries  in  the  Table  refer  to  the  following 
posi cions: 

PI  :  Wing  Tip;  Top 
PI*;  Wing  Tip;  Bottom 
P2  :  Wing  Center;  Top 
P2‘:  Wing  Center;  Bottom 
P3  :  Wing  Root;  Top 
P4  ;  Fuselage 

For  simulating  the  aircraft  over  ground  plane  configuration,  the 
model  was  placed  an  eguivalent  distance  from  the  ground  plane  as  if 
it  were  sitting  on  its  wheels.  This  placed  the  bottom  of  the 
fuselage  at  the  wings  2  cm  from  the  ground  plane. 

7.0  Comparison  with  University  of  Michigan  CW  Data  and  with 
Analytical  Predictions 

The  University  of  Michigan  conducted  CW  tests  of  three  test  objects 
[Ref.  5],  Two  of  these  tests  (i.e.,  cylinders  and  a  scale  model  of  a  747 
aircraft)  correspond  to  tests  described  in  this  report.  In  this  section, 
the  two  sets  of  data  will  be  compared.  Enough  detail  of  the  Michigan 
set-up  and  data  will  be  given  here  to  make  the  comoarisons  meaningful. 
Interested  readers  should  consult  the  original  Michigan  report  for 
additional  information. 

The  Michigan  measurements  include  induced  current  and  charge  results  for 
plane  electromagnetic  wave  incidence.  These  were  performed  in  the  model 
frequency  range  of  45  MHz  to  4.25  GHz  using  the  image  plane  technigues. 


5.  Liepa,  V.  V.,  et  al.,  "Surface  Field  Measurements  with  Imaqe  and  Ground 
Planes",  AFWL  Sensor  and  Simulation  Note  244,  Albuquerque,  NM,  November, 
1977. 
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TABLE 

4.  SUMMARY 

OF  TIME 

DOMAIN  RESULTS 

FOR  747 

SCALE-MODEL 

AI PCPAFT 

Aircraft 

Location 

F  ree 

Space  (FS) 
or 

Ground 

Plane  (GP) 

D-Sensor 

Fig.# 

r /r inc 

F  i  q .  * 

• 

B- 

Peak 

Sensor 
'iq.#  - 

H/H  1  r>r- 

iq.  8 

°eak 

PI 

FS 

FI 

F5 

13 

N/A 

N/A 

N/A 

PI  ‘ 

FS 

F 1 3 

F  1  7 

1? 

N/A 

N/A 

N/A 

GP 

FI  1 

ns 

90 

N/A 

N '  A 

N/A 

P2 

FS 

F2 

F  7 

7 

F3 

rq 

P2 ' 

FS 

N/A 

N/A 

N/A 

F  1 4 

"19 

GP 

N/A 

N/A 

N/A 

F12 

r  16 

2  S 

P3 

FS 

N/A 

N/A 

N/A 

F  4 

c9 

S 

P4 

FS 

N/A 

N/A 

N/A 

FS 

HO 

4 

N/A  -  Data  not  taken. 
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Data  were  obtained  for  free  space  and  in  the  presence  of  perfectly 
conducting  ground  planes.  The  only  Michigan  data  quoted  here  correspond 
to  cases  where  a  one-to-one  comparison  could  be  madp  with  LLNL  data. 

Comparisons  between  the  cylinder  data  will  be  taken  up  first  and  the 
results  summarized  in  tabular  form.  The  Michigan  and  LLNL  cylinder  radii 
and  length  are  very  close  (radii;  5.08  cm  vs.  5.05  cm  lengths;  50.83  cm 
vs.  50.5  cm).  The  positioning  of  the  sensors  is  not  quite  close  (0  sensor: 
4.83  cm  vs.  10.1  cm  from  top  plate;  8  sensor:  4.83  cm  vs.  2.5  cm  from 
bottom  plate;  refer  to  Figure  6).  It  is  not  obvious  how  much  uncertainty 
this  introduces  into  the  comparisons.  Most  probably  the  effect  on 
surface  current  measurements  near  the  bottom  of  the  cylinder  may  be 
small;  the  uncertainty  may  be  larger  for  the  D  measurement  close  to  the 
top  of  the  cylinder,  for  the  cylinder  over  qround  olane  measurements, 
the  height  h  in  all  cases  are  almost  identical  (10  •  cm  vs.  10.16  cm; 

50  •  cm  vs.  50.8  cm;  100  cm  vs.  101.6  cm;  and  200  cm  vs.  203.2  cm).  The 
results  of  these  comparisons  are  shown  in  Table  5  where  the  Deak  value  of 
the  transfer  function  at  the  fundamental  resonance  is  tabulated. 

For  the  747  scale-model  aircraft,  the  comparisons  are  summarized  in 
Table  6.  For  details  on  location  notation,  please  refer  to  Section  6.3. 

In  addition  to  Michigan's  CW  data  for  the  cylinder,  analytical 
predictions  were  made  for  the  same  configuration  by  Sancer,  et  al. 

The  predictions  are  entered  in  Table  5  in  the  third  column.  Ref.  6  also 
included  predictions  for  h  =  7.5  cm  and  25  cm  where  transfer  function 
peaks  at  the  fundamental  frequency  of  404  and  67,  respectively,  are 
predicted.  In  order  to  facilitate  the  comparisons,  the  LLNL,  U.  of  M. 
and  the  Ref.  6  predictions  for  the  peak  transfer  function  at  the 
fundamental  frequency  are  shown  in  cigure  14  for  the  cylinder  above 
ground. 

In  Fiqures  15,  16,  and  17,  actual  frequency  domain  spectral  content  are 
shown  for  various  cylinder  conf igurations.  These  include  LLNL  and 
University  of  Michigan  measurement  predictions  and  analytical  predictions 
in  Ref.  6.  The  top  curve  in  each  figure  includes  the  smoothed  out 
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Table  5. 

Summary  of  the  Comparisons  for  the  Cylinder 
(Transfer  function  peak  value) 

Transfer 

Function 

Height 

h 

Annie 

9 

U.  of  M.  LLNL 

Analytical  Prediction 
(Ref.  6) 

TABLE  6.  SUMMARY  OF  COMPARISONS  FOR  THE  747  SCALE-MODEL 
AIRCRAFT  (TRANSFER  FUNCTION  PEAK  VALUE) 


Location 

Free  Space 
(FS) 
or 

Ground  Plane 
(GP) 

Transfer 

F  unct i on 

T  ype 

IJ.  of  M. 

L  .  N . 

PI 

FS 

E/F inc 

16.0 

13 

P2 

FS 

H/n'hc 

9.6 

■> 

GP 

H/H 'nC 

16.7 

76 

P3 

FS 

H /H 1 nC 

7.0 

P4 

FS 

H/H>"C 

8.0 

4 

PI  -  Winq  Tin,  Too 
?2  -  Winq  Center,  Top 
P3  -  Winq  Root,  Top 
P4  -  Fuselaqe 


version  of  the  Michiqan  data  as  reported  in  Ref.  f>,  ind  the  analvfir->' 
Dredictions.  The  bottom  portion  is  the  correspond  1  no  LLNL  data.  rne 
free  space  case  is  shown  in  Fiqure  IS.  The  case  for  a  cylinder  ’6,4  --i 

over  a  perfectly  conductinq  qround  o 1 ane  is  shown  in  Fiqure  16.  Toe 

results  for  a  cylinder  very  close  to  oround  are  shown  in  figure  !?.  'h 

Michiqan  data  and  the  predictions  are  for  h  -  7.6  ■-«),  The  corresp on.:  i n 

LLNL  data  is  for  h  =  8.?  cm  since  the  cylinder  could  not  be  brought  ,nv 
closer  to  the  qround  plane  because  of  nhysical  obstructions  at  the  n-:se 

In  Fiqure  14,  note  that  for  separations  larqpr  than  u  -  oO  c.m,  both  too 
LLNL  and  University  of  Michiqan  data  are  cIosp  and  exhibit  similar 
trends.  For  separations  closer  than  h  =  60  cm,  the  LLNL  data  exoiM*, 
larqer  transfer  function  peaks  than  the  correspond inq  University  e‘ 


6.  Sancer,  M.  I.  et  al.,  "Formulat  ion  of  r 1 pctromaqnet  n;  nulse  r*ter.i>’ 
Interaction  Above  a  Lossy  carth7Comnarisnn  of  Numerical  Results  w'*n 
Experimental  Data  for  Limit! nq  Cases,"  AF’WL  Interaction  Notes,  Note  ',4 
Albuquerque,  NM,  October  1978. 


T'tichiqan  data.  In  this,  the  LLNL  data  exhibit  the  trend  suggested  by  the 
analytical  oredictions  for  h  <  10  cm.  There  is  a  wide  discrepancy. 

however,  at  h  =  25  cm.  Considering  the  difficulties  in  takinq  such 
hiqhly  resonant  data,  this  may  not  be  too  surorisinq.  The  aqreement 
between  the  LLNL  and  the  University  of  N'ichiqan  data  for  the  747 
scale-model  aircraft  is  also  quite  encouraging  as  shown  in  Table  5. 

Tne  spurious  neeks  observed  at  the  hiqher  frequencies  in  the  LLNl  dat i  in 
figures  15  and  17,  which  do  not  annear  in  the  University  of  Micninan  fata 
are  produced  by  the  division  of  one  transform  by  another.  The  small 
baseline  numbers  generated  by  the  tYT  routine  are  very  noisy.  When  •ovgo 
small  baseline  numbers  approach  zero  the  division  process  will  cause  > 
numerical  noise  neak  to  appear  since  J/X  goes  to  a  large  value  as  < 
approaches  yprp. 

S.O  r*o  1  e  extraction  Results 


Three  parameter  estimation  algorithms  were  aopliod  to  the  data  to  t 

the  poles.  These  are  the  Prony  f echn i ques , *  extended  <alman  filler,  ’-r 
and  the  maximum  likelihood  estimator.  All  three  are  discussed  <vm 
compared  in  terms  of  their  respective  theories  in  Append  i  *  "  .  r,m 
application  of  the  rKF  and  the  maximum  likelihood  estimator  t'  *  o.  <  .*  ■ 
are  discussed  in  Apnendi*  and  the  results  tabulated  in  this  set '  ’  . 

Tne'e  two  p,ir  amer  er  est  n-ators  w<*r«-  •an  1  y  ippl  md  t ..  the  ry’w  !•■  . 

'he  jppl  icat  ion  of  the  nrony  te  hniqi"  n  revere  1  m  knoo>.<'v  T.  ' 
algorithm  war  applied  to  all  three  objects  and  .  ’  I  f*i.<  results  ?r.» 
tabulated  in  that  appendix.  -  •*  Purposes  nf  o-’u-ir  i  •.•••>  wi*h  t  ••  <  ■  •• 

t  furore  t  ical  ;>redi  lion-,  and  I'.viip  f  roa  'no  » .i»n,»r  e\*  •  .♦  • 

*’  r  on  y  results  for  tne  cylinder  p..  extracted  ( r  v  *n.-  amvn  1 1 . 
in  tabulated  fashion  in  this  i  m. 


Tne  ;,rony  technique*  could  really  '«■  •■,r-'ed  '•  .  uter  a  '  ive"  *r  my. 
many  runs  through  'he  data  were  \  imp  I  i  shed ,  i!l  .  1  ' 

niise  and  i n<  rease  si  dial  levels. 


i  ) 


*  i  M  .  * 


Figure  14. 


Predictions  tor  the  cylinde 


The  result  of  the  oole  extraction  exercise  is  shown  in  Table  7.  Here, 
the  Prony,  Extended  Kalman  Filter,  and  Maximum  Likelihood  Estimator 
oredictions  are  listed  as  well  as  analytical  nredictions  by  Shumpert  and 
Galloway  [Ref.  7].  Only  the  cylinder  case  is  considered  for  three 
configurations:  free  space,  h  =  50  cm  and  h  =  10  cm.  Only  the  first 
three  poles  are  shown.  The  table  shows  estimates  for  the  poles  (real 
part,  cr,  and  imaginary  part,  oj)  .  The  maximum  likelihood  identifier  was 
run  using  two  different  underlying  signal  model  assumptions  tyoes  A  and 
8.  These  two  model  types  are  decribed  briefly  in  the  text  below. 

We  make  the  following  comments  concerning  the  tabulated  results.  The 
second  harmonic  (pole  #2)  parameters  are  generally  nard  to  identify 
because  very  little  second  harmonic  information  is  in  the  measurement 
data.  The  current  probe  on  the  cylinder  happened  to  be  near  a  node  for 
the  even  harmonics.  Therefore,  the  Prony  technique  did  not  even  assign  a 
pole  at  the  second  harmonic.  The  maximum  likelihood  method  forces 
identification  of  second  harmonic  parameters  because  the  system  model 
internal  to  this  algorithm  requires  that  the  second  pole  be  present.  The 
high  estimated  dampinq  terms  (<r)  can  be  explained  when  we  consider  that 
second  harmonic  exists  in  the  initial  transient  (the  electromagnetic 
impulse)  but  is  not  above  the  noise  level  in  the  remainder  of  the  data 
record.  Analytical  predictions  for  the  third  harmonic,  50  cm  and  10  cm 
data,  are  not  available. 

Because  second  harmonic  information  is  virtually  unobservable  in  the 
data,  we  restrict  our  comparisons  to  the  fundamental  and  third  harmonic 
only.  Fundamental  frequency  parameter  estimates  usinq  analytical,  Prony, 
and  maximum  likelihood  (model  tyoes  A  and  B)  methods  are  plotted  in  the 
s-plane  in  Fiqure  18.  The  graph  shows  the  migration  of  the  pole  (the 
fundamental)  as  the  distance  between  tne  cylinder  and  the  wall  becomes 
smaller  according  to  the  results  of  each  analysis  method.  Fiqure  19 
shows  a  similar  graph  for  tne  third  harmonic  pole. 


Shumpert,  T.  H.,  and  Galloway,  0.  J.,  "Transient  Analysis  of  a  Finite 
Length  Cylindrical  Scatterer  Very  Near  a  Perfectly  Conducting  Ground" 
AFWL  Sensor  and  Simulation  Notes,  Note  226,  Albuquerque,  New  Mexico, 
August  1976. 


38 


1 


TABLE 

7.  SUMMARY 

OF  POLE 

LOCATIONS  FOR 

CYLINDER  (s  = 

O  +  j/.;) 

Pole  4 

Ana  lyt ica 1 
Prediction 

Prony 

Extended 

Ka Iman 

F i 1  ter 

Max . 

Likelihood  A* 

Max . 

Likelihood  3** 

Free  Space 

- 131.4 

-99.6 

-190 

-215 

-137 

775.4 

720.7 

750 

526 

829 

1  h  =  50  cm 

-66 . 3 

-74 

-65 

-102 

-70 

736.2 

541.0 

736 

7  70 

722 

h  =  10  cm 

-16.9 

-22.9 

-25 

-28.9 

-25 

324.0 

661 .0 

758 

767 

775 

Free  Space 

-206.4 

-N/A- 

-N/A- 

-654 

-N/A- 

1665.3 

-N/A- 

1052 

2  h  =  50  cm 

-301.5 

-N/A- 

-N/A- 

-1536 

-N/A- 

1539.0 

-N/A- 

1540 

h  =  10  cm 

-55.3 

-N/A- 

-N/A- 

-723 

-N/A- 

1653.0 

-N/A- 

1534 

Free  Sp^e 

-236 

-297.3 

-710 

-721 

-654 

2523.6 

2433.5 

2250 

1578 

2480 

3  h  =  50  cm 

-N/A- 

-258.9 

-460 

-440 

-206 

-N/A- 

2333.2 

2208 

2310 

2165 

h  =  10  cm 

-N/A- 

-53.7 

-160 

-50.5 

-200 

2311.6 

2274 

2301 

2379 

Note: 

For 

*  A: 

W  i  th 

**  B: 

W  i  th 

N/A: 

data 

bottom  oj. 

noise  models  included. 
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The  FKF  and  maximum  likelihood  B  results  appear  to  best  agree  with  the 
analytically  predicted  values. 


1 


A  significant  difference  between  two  types  of  internal  signal  models  (A  ^ 
B)  used  should  be  pointed  out  at  this  time.  As  applied  to  the  pole 
extraction  problem,  the  maximum  likelihood  A  internal  system  model 
assumes  that  the  cylinder  has  been  illuminated  by  an  ideal  impulse,  and 
the  resulting  response  is  thereby  a  sum  of  damped  exponetials.  The  EkE 
and  the  maximum  likelihood  R  models  have  the  more  qeneral  assumption  that 
the  cylinder  is,  over  the  period  of  the  experiment,  randomly  excited. 

The  maximum  likelihood  A  identifier  therefore  attempts  to  fit  the  data  to 
a  sum  of  exponentials,  whereas  the  EKF  and  maximum  likelihood  R 
essentially  fits  the  autocorrelation  function  (ACT)  of  the  data  to  tne 
ACF  of  a  sum  of  expontials. 

It  is  not  yet  clear  which  technique  is  best  for  use  on  the  EMD  problem, 
since  EM?  data  is  usually  a  combination  of  responses  to  an  ideal  impulse 
and  to  random  excitations.  A  future  task  would  be  to  include  botn  the 
ideal  impulse  and  random  excitation  models  into  EKF  and  maximum 
likelihood  identifiers.  This  would  fold  in  to  the  proposed  statistical 
characterization  of  the  EM?  range  as  described  in  Appendix  C. 

9.0  Conclusions  and  Recommendations 


In  this  work,  we  have  demonstrated  the  utility  of  the  Transient  Ranqe 
Simulator  to  record  and  process  transient  data.  Three  objects  were 
tested,  namely  a  cylinder,  crossed  cylinders,  and  a  scale-model  747 
aircraft.  The  data  were  processed  to  find  the  frequency  spectral  content 
and  the  poles.  In  each  case,  comparisons  were  made  with  other  available 
test  data  or  analytical  predictions.  The  following  paragraphs  contain 
pertinent  comments  concerning  several  aspects  of  this  work. 

A  recurring  problem  in  time-domain  measurements  is  the  lack  of  stable, 
hiqh-ampl itude  pulse  generators.  This  series  of  experiments  used  the 
I KOR  IMP  impulse  generator,  which  produces  an  impulse  of  voltage,  300  os 
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Analytical  prediction 
Interactive  Prony 
Extended  Kalman  filter 
Maximum  liklihood  A 
Maximum  likelihood  B 


Figure  18.  Migration  of  the  Fundamental  Pole. 
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•  Analytical  prediction 

- • -  Interactive  Prony 

.  Extended  Kalman  filter 

—  - -  Maximum  liklihood  A 

- Maximum  likelihood  B 


Figure  19.  Migration  of  the  third  Harmonic. 


wide  at  6(]%  amplitude  point  into  a  50-f2  load  (Ref.  2).  Other 
commercially  available  pulsers  do  not  nave  the  amplitude  or  soectral 
content  required  to  make  wide-band  measurements.  Coupled  to  this  problem 
are  the  sensors  used  to  make  the  measurements.  Ideal  sensors  have  nigh 
sensitivity  and  are  sufficiently  .mall  in  size  to  minimize  the 

perturbation  of  the  fields  they  are  designed  to  measure.  cor  these 

•  • 

measurements  we  used  the  AFWL  ACD-1  0  and  the  M6L-8  h  sensors  -  tw 
wel 1 -desiqned  devices,  flat,  even  with  the  highest  amplitude  pulsar 
available,  we  found  ourselves  close  to  the  noise  level  on  several 
measurements.  Advances  will  probably  have  to  come  in  the  design  of 
pulsers  if  measurements  like  those  employed  here  are  to  be  extended  *-n 
more  stringent  conditions. 

The  possible  use  of  a  hiaher  amplitude  pulse  venerator  to  improve  mr 
signal  to  noise  ratio  is  linked  to  the  question  of  signal  process’' 
methods.  If  higher  amplitude  signal  generators  with  the  proper  f r^g  r*n-. , 
spectrum  cannot  be  luxated,  then  improved  signal  processing  met ho  Is  nee  : 
to  be  used  to  compensate  for  the  ooor  signal  to  noise  ratio. 

The  greatest  difficulty  encountered  in  this  project  was  the  recoriin  :  in 
processing  of  the  highly  resonant  responses  of  the  objects  as  thev  were 
brought  closer  to  the  ground  plant.  This  necess i tai.e<l  t  < i n q  longer  time 
records  which  accentuated  problems  of  sampling,  drift  and  reflective. 
These  in  turn  cause  problems  in  the  determination  of  the  spectral  content 
by  classical  means  such  as  the  Fast  'ourier  Transform  or  the  F i Ion 
Metnod.  Most  of  these  problems  were  solved  by  recording  the  data  with  a 
Tektronix  791?  in  conjunction  with  an  LSI-11. 

One  of  the  recurring  problems  durinq  the  course  of  this  project  >*,.is  th" 
lack  of  adequate  character i zat ion  of  the  range,  a  successful  anpl irat ion 
of  signal  process inq/ident i f icat ion  methods  to  the  range  will  require  one 
to  thoroughly  characterize  the  nature  of  the  measurement  system  and 
background  noise.  Current  improvements  on  the  range  instrumentation 
should  improve  the  attainable  siqna l-tu-noise  ratios  significantly; 
however,  specific  tailoring  of  the  measurement  techniques  to  signal 
processing  need,  will  need  to  tie  implemented  to  assure;  II  enough  oati 
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are  taxen  to  assure  observability  of  the  parameters  of  interest,  and  2) 
data  rates  are  sufficient  to  prevent  aliasing.  A  further  subject  of 
interest  is  the  characterization  of  an  HMP  resoonse  by  a  more  complex 
dynamic  model  than  the  simple  sum  of  damped  sinusoids  model  which  is  now 
used.  This  model  would  include  generally  nonlinear  coupled  differential 
equations  and  correlated  driving  noise.  The  complex  model  would 
hopefully  incorporate  knowledge  of  the  phenomenological  asoects  of  the 
experiment  such  as  elements  of  the  EM  field  equations,  pulse  reflection 
off  the  ground  plane,  etc. 

In  analyzing  the  data  and  obtaining  the  transfer  functions  for  the 
various  objects,  FFT  techniques  were  used.  In  retrosoect,  this  seems 
like  a  rather  Door  method  to  perform  the  signal  processing.  The  process 
of  dividing  frequency  transforms  is  prone  to  large  errors.  We  recommend 
that  in  the  future,  time  domain  signal  processing  methods  which  do  not 
require  EFT'S  be  used.  Two  such  methods  are  becoming  currently  available 
and  have  been  discussed  in  the  main  body  of  the  report;  namely  the 
extended  Kalman  filter  and  the  maximum  likelihood  estimator.  These 
process  the  time  domain  data  directly.  Models  of  the  incident  pulse  can 
be  employed  to  "remove"  its  effect  from  the  data.  In  this  way,  the 
numerical  problems  created  by  frequency  domain  (FFT)  division  are  avoided. 
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Append  l  <  •> 

i.;.Nl  Frans  ient  tange  F  ar  i  1  i  f  y 

Measurements  at  the  tLNL  Transient  -Tange  -a  ility  are  made  over  i  la-'q- 
x  8.5  m  aluminum  covered  ground  olane  ■loun t •  )  J  feet  •.iff  tne  floor  using  i 
conical  antenna  whicn  extends  from  tne  image  olane  to  the  •"•ooi',  oilinq,  1  i 
instrumentation  is  located  beneath  the  image  olane,  and  at  this  facilit/  a 
sampling  oscilloscooe  is  utilized  as  the  siunal  sampling  device. 

The  operation  of  the  transient  range  is  conceptually  very  simple.  A  very 
narrow  electrical  pulse  is  used  to  drive  a  wide-bandwidtn  radiating  ant-no  i 
which  illuminates  tne  target  under  test.  The  surface  densities  of  currents 
and  charges  induced  on  the  target  are  measured  by  surface-mounted  i  an  :  > 
sensors.  The  outputs  from  the  sensors  are  fed  to  a  sampling  osc i 1  lose om’ 
which  records  the  response  of  the  taroot  as  a  function  of  time.  The  rr.su  I  tr- 
may  be  disolayed  in  the  time  domain  or  transformed  to  the  frequency  do-"  a  in  via 
the  Fourier  transform.  Figure  A1  shows  a  block  diagram  nf  the  system;  the 
diagram  shows  that  a  minicomputer  controls  tne  ooeration  of  the  equipment  and 
handles  the  data  loqqing  tasks.  The  nature  of  the  samplinq  oscilloscope  is 
such  that  not  one,  but  many  pulses  must  be  radiated  by  the  source  antenna  to 
provide  a  complete  time-history  waveform.  The  computer  controls  the 
oscilloscope  and,  in  the  orocess,  obtains  512  equally  spaced  samples  of  tne 
transient  waveform.  These  transient  data  are  then  plotted  on  an  on-line 
plotter  and  then  written  onto  magnetic  tape  for  later  off-line  data  processing. 

Two  factors  determine  the  working  oandwidth  )f  the  transient  measurement 
system.  At  tne  low  frequency  end,  the  rang  clear  time,  or  time  it  takes 
unwanted  reflections  from  discontinuities  in  the  cone,  and  from  the  walls  and 
surrounding  environment  to  reach  the  target  area,  establishes  a  lower  limit. 

For  the  LLNl  ranqe,  reflections  from  ,imularor  edges  (building  ceiling  and 
other  objects  in  the  room  have  been  found  to  be  nondetectaDl e)  and  the  eanqe 
have  been  used  to  obtain  data  records  in  excess  of  200  ns.  For  limited  number 
of  samples,  this  limits  the  time  •'•<»solut ion  and  may  introduce  large  sampling 
errors.  A  200  ns  time  record  corresponds  to  a  minimum  useful  frequency  of 
5  MHz.  At  the  high  frequency  end,  the  limits  are  a  complex  combination  of 
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Figure  AT.  Schematic  of  LLNL  Transient  EM-Measurement  Facility. 
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with  the  do  1  ser  output  as  a  t  r  t  ■  i  • )  •  *  r  ,  .  i  r  ■ .  ’nr  *  ne  s  /s  t  ."r .  ",j  . 
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1.  Asxs  the  experimenter  for  v  .  ■»  ),»•..  ■- 1 r  i  ,jn  s  an.)  scope  - 1 
f actors . 

\  Jlots  results  l  indue  inn  s:  -l.»s  ami  ritUgi  or.i,, 

plotter. 

3.  Records  acfiiinu  1  ated  experi  m-rit  a  i  i  )■  i  on  magnet  :  »  <  * 

oft -line  post  processing. 

4.  Controls  the  sampling  use  l  1  I  os -,>;■>  mi  .-ruing  n 
and  averaging  of  t fie  sample.:  .leim  ,. 

fne  above  data  collection  system  was  :  murid  r- , ,  g,-  adequate  •  m  t  in. 
in  excess  of  ?i)  ns.  oir  longer  records,  **«•  •> /  > t **im  was  ..•*  .%  i  *  ■ 
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to  control  and  Drocess  tne  data  was  t  .and  t  >  lie  far  ,.i;i..ri 
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several  t  y  p  •  *  s  v  .1 1  unfed .  ■  iqure  A?  shovs  the  transfer  functions  of  the  lin<-s 

tested.  Of  these,  the  7/8  in.  Hr  I i ax  Tektronix  unit  performed  best  and  was 
tnerefore  selected  for  tnese  measurement. 

Conical  Antenna--The  final  item  we  improved  was  the  conical  antenna  used  to 
radiate  the  transient  CM  ou'se.  By  increasing  the  cone  anqle,  as  shown  in 
Figure  A3  [Ref.  Al],  the  conical  antenna  can  be  marie  to  match  the  S0-J2 
feed-line.  This  occurs  at  a  cone  half-angle  of  =:  43°. 

To  evaluate  th  desinn  of  the  cone,  we  first  utilized  a  r j re-Oomain 
Ref lectomete--  VTPR)  to  match  the  cone  impedance  <-n  toe  fnedline.  -inure  "4 
shows  a  cross  section  of  the  transition  region  and  the  wav  the  ted  is 
adjusted  for  an  impedance  match. 

By  adjusting  the  length  of  the  feed  screw,  we  can  raise  and  lower  the  rone 
into  the  cylindrical  coax  section,  [f  me  cone  tip  is  too  long,  the  match  to 
the  cone  is  inductive;  if  it  is  too  short,  th.-  capacitance  of  the  cone  to 
sides  of  the  coaxial  feed  can  be  observ'd.  nrooer  adjustment  of  the  feed 
length  can  be  attained  for  a  perfect  match.  This  is  shown  on  the  TPR  plots  of 
F igure  A5. 

Away  from  the  feed  region,  the  cone  is  composed  of  three  reqions,  as  shown  in 
Figure  A6.  The  bottom  section  is  formed  with  a  niece  of  sheet  brass,  while 
the  upper  sections  are  comoosed  of  16  symmetrical  1 v  located  wires  whicn  extend 
from  tne  brass  cone  to  the  room  ceiling.  Aluminum  foil  was  also  used  to  cover 
the  lower_part  of  the  wire  section.  Above  the  foil,  each  of  the  18  wires  are 
loaded  with  eight  discrete  lrJ0-fi  resistors.  These  leads  help  to  attenuate  the 
pulse  propagating  on  the  cone,  thus  minimizing  the  pulse  rpfierted  from  the 
ceilinq  back  into  the  experimental  area.  Fiqure  A7  shows  a  TPR  plot,  pf  the 
whole  conical  structure  in  which  the  various  sections  of  the  cone  may  he  noted. 


Al.  Antenna  Engineering  Handbook,  Jasik,  H.,  r<l. ,  McGraw-Hill  Book  To., 
New  York,  1961,  pp.  3-11. 
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Figure  A3.  Cnaracterr.t  i>  I'npegjruie  o'"  .ironical  -1  i  ;>o  1  r.  v 
cone  angl e  [Ref .  ?) . 
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Figure  A4 .  Cross  section  of  Monocone  feed  region  of  the  conical 
antenna  showing  the  adjustment  screw  used  to  raise 
and  lower  cone  in  the  coaxial  fitting. 


iqurp  A 5.  T'iR  of  «on*»  Tcq ion  showinq  prooor  .id.i : -M'r.t  of  ;‘v’p'i  1  n a t » 


Figure  A6.  Construction  details  of  monocone  antenna 


Figure  A7.  TDR  Plot  of  monocone  input  impedance,  showing 
the  various  regions  of  the  cone. 


Append  i  x  'I 

Rationale  For  Cylinder  Longtn  Choice 


Several  factors  influence  the  choice  of  a  te*t  object  size  for  scale  model 
measurement,  and  the  selection  is  a  trade-ofr  between  these  various  factors 
discussed  in  the  following  paragraphs. 


Phase  Error;  The  incident  wave  that  illuminates  the  test  object,  at  too 
LLNL  Transient  Range  facility  is  a  spherical  wave.  Therefore,  the  whole 
object  is  not  illuminated  at  the  same  tone.  This  gives  rise  to  pnase 
errors  that  depend  on  the  frequency  components  of  interest. 


For  a  test  object  that  is  a  vertical  cylinder  of  length  ,/2  (see  "ion-. 
5)  located  a  distance  R  from  the  base  or  the  monocone  antenna,  the 
maximum  time  difference  over  the  object  in  the  arrival  time  of  the 
incident  pulse  is 


At 


+  L2/4  -  R  1 


which  for  R  >>  l/2  reduces  to 

At  =  L “ /8  Rc 


where  c  is  the  speed  of  light.  At  the  fundamental  resonance  of  the  test 
object  ( X.  =  21),  this  transforms  into  phase  error  A0,  where 


+  L2/4  -  R  1 


(<33) 


or  A0=  n  L/8  R.  [R  L/2]  .  (34) 

For  the  cylinder  used  in  toe  experiments  reported  here,  i_  =  1.01  m  and 
R  =  2.43  m,  which  gives  a  A0~9°.  Tiv  statement  of  work  called  for  a 
phase  error  not  to  exceed  45°  at  the  f  .  idamental  resonance.  Mote  that  if 
I.  is  made  to  increase,  the  phase  ’rror  vill  similarly  increase  almost 
linearly  with  L.  Therefore,  it  i,  ore' arable  to  keen  the  larqest 
dimension  of  the  test  obiect  as  small  . oossiole. 


Oar  vr- 


Amplitude  Variation;  Another  reason  for  keeping  L  small  is  to  avoid 
large  incident  field  amplitude  variations  across  the  object.  In  ratio 
this  is  given  by 

Amplitude  Variation  =  +■  L^/4  -  R  ]  /R  (35) 

which  for  R  >>  l/2  becomes 

Amplitude  Variation  =  L^/8  8^  (36) 

For  the  cylinder  in  question,  this  is  ^  2.%. 

Signal-to-noise;  This  is  a  function  of  the  pulser  energy  and  bandwidth, 
target  response  and  sensor  sensitivity.  If  everything  else  is  left 
uncnanged,  better  signal-to-noise  ratio  is  usually  obtained  by  larger 
sized  objects.  A  larger  sized  object  is  also  easier  to  instrument.  As 
we  have  seen  however,  this  can  only  be  pushed  so  far,  since  phase  and 
amplitude  errors  will  become  unacceptable. 

It  was  decided  at  the  start  of  the  effort  that  a  good  compromise  for 
object  size  was  a  largest  dimension  equal  to  about  .5  m. 


Appendix  C 

Pole  Estimation  Techniques 

1 .0  Introduction 

Tne  study  of  E^P  phenomena  has  promoted  tin*  develjpment  of  techniques  t > 
investigate  transient  electromagnetic  response  data.  The  characterization 
of  EMP  transient  response  information  is  a  matter  of  national  concern. 
Since  large  amounts  of  data  are  necessary  to  nointwise  define  an 
arbitrary  transient  response,  it  is  quite  reasonable  to  "identify"  i 
parameteri zation  or  model  of  the  "response".  Tne  model  developed  is 
useful,  not  only  to  merely  parameterize  the  response.  Put  also  to  nive 
more  meaningful  information  about  the  pnysical  process  producing  tne 
response  itself. 

In  this  appendix  we  discuss  the  implementation  of  some  siqnal  processing 
algorithms  which  can  be  used  to  "estimate"  the  parameters  of  an  electro¬ 
magnetic  response  model  from  noisy  transient  measurements.  Tne 
techniques  employed  range  from  simplified  algoritnms  wnich  perform  well 
for  hiqh  siqnal-to-noise  ratios,  to  complex  model -based  estimator-;,  wh-cn 
perform  well  for  low  siqnal-to-noise  ratios.  In  Section  ?  we  present  the 
necessary  background  i rif orination.  The  various  algoritnms  are  discusse' 
(simply)  in  Section  3.  In  Section  4,  tne  application  to  transient 
data  is  presented. 

2.0  Sacxground 

In  electromagnetic  wave  theory  it  is  imssible  to  represent  the  rev>->n-,- 
of  an  object  to  various  excitations  ;>y  tne  singular!!  /  expansion  ’’“'ho; 

( St£M ) .  The  SEM  represents  an  electro-i-a  met  ic  variable  'fieiq,  curro-it., 
etc.)  as  tne  impulse  response  of  the  obi  eel  (_3ef.  11  j,  i.e., 

•> .  t 

2-  h^e^pv^r)-  '  '  1  ' 


h  / 


Vi*1’  = 


whore 


Up  vector  Vpjlse  response 

complex  coupling  coefficient 

v^  complex  natural  mode  describing  the  behavior  of  Up  over  the 
object 

e  exciting  field  character i st ics  (e.g.,  pol ar i zat i on,  direction 
of  evidence,  etc.) 

r  spatial  coordinates  or  position 

s-  complex  natural  frequency  (or  oole,  or  natural  resonance) 

The  sets  of  parameters  ({s^},  { _v - ( jr) })  are  dependent  on  the  object 
parameters  only  and  independent  of  the  excitation.  The  effect  of  the 
exciting  wave  is  contained  entirely  within  the  set  of  couplinq  coefficients 
{■e,  s.)  which  are  independent  of  the  position  on  the  body.  Thus,  the 
electromagnetic  interaction  is  completely  characterized  by  these  sets.  In 
fact,  the  response  to  an  arbitrarily  shaped  waveform  can  be  generated  usinq 
concepts  of  linear  system  theory  where  the  response  y  to  an  arbitrary 
impulse  is  given  by  the  convolution^) 

y(j\t)  =  Up(r,t)  *  u(r,t)  (33) 


Implicit  in  (C2)  is  that  e  is  the  same  for  the  new  excitinq  waveform. 
However,  ({s^,  {v^(r)})  are  invariant;  therefore,  these  sets  can  be  used 
to  "parameterize"  a  given  object  for  any  excitation.  We  are  not  concerned 
at  this  ooint  in  the  partitioninq  of  the  natural  modes  and  couolinq  coeffi¬ 
cients,  so  we  define  the  set  of  complex  residues  at  a  point  r  as 


C-jdo):  =ni  (e,si)v.(ro) 
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and  for  this  work  concern  ourselves  only  with  scalar  response  functions.'*’ 
Thus,  the  impulse  response  of  the  linpa>"  system  of  (~°)  can  be  represented  as 

^  It  should  be  noted  that  the  sophisticated  node! -based  estimators  discussei 
subsequently  can  be  used  to  identify  separately  the  .  and  v.  parameters, 
if  desired,  as  well  as  vector  response  functions  (multiple  measurement, 
instruments);  however,  this  work  was  not  feasible  in  the  il  lotted  time. 

vl 


where 


y(t)  =  h  ( t ) 


N 

*  6  ( t)  =  Z 
i  =  l 


C  .e 

1 


V 


(CM 


s.j  ;  =  <r-  +  j  w-  ,  <r  the  damn  ion  ratio  and  u.  the  natural  frequency 

H( t )  is  the  object  impulse  response  at  position  r^;  >.p.t  :|  (r^.t). 

Thus,  the  parameter  i zat  ion  of  th-’  ob  ject  can  he  stated  simply  by  the 
electromagnetic  parameter  estimation  problem. 

"Given  a  set  of  noisy  electromagnetic  response^  measurements  z{ t)  , 
find  the  (best)  (minimum  variance)  estimates  ({a.,  w.},  (c-j)  character¬ 
izing  an  unknown  object." 

We  will  assume  that  the  noise  contaminates  the  response  y  as 

z(t)  =  y(t)  +  v(t)  f 

where  v  is  zero  mean,  Gaussian  with  covariance  1. 

Refore  we  begin  discussina  the  various  estimation  algorithms  aopl’^d  t 

the  problem,  we  must  definp  an  alternate  way  nf  representing  a  linear 

system  which  is  equivalent  to  (CD  and  (C4).  decal  1  from  ordinary 

differential  equations  ’Ref.  C?]  that  ('41  represents  the  solution  nf  a 
f  h 

M'  order  differential  equation.  It  is  well  known  that  this  equation 
can  be  broken  down  to  the  solution  of  ')  first  irder  differential 
equations  of  the  general  form; 

D  t )  -  F  x_(  t)  e  q  u(  t) 

y(t)  -  hrx(t)  DM 

*  This  representation  is  not  limited  only  t  .  a  I,*.-  >.  vs  t  «*•  > ;  e .  q . ,  u  am1  y  ;an 

r 

be  vectors  and  q,  h1  become  matrices. 


where 


x  is  the  N-state  vector,  u,  y  are  the  respective  input  and  output. 

F  is  a  N  x  N  matrix  and  q,  h  are  N-vectors 

This  representat ion  is  called  the  "state  space"  form  in  linear  system 
theory  [Ref.  C2]  and  forms  the  basis  of  various  parameter  estimation 
schemes  [e.g.,  see  Ref.  C4].  It  is  easily  shown  that  the  impulse 
response  of  (C6)  is 

r  ft  ^  5  i 

y(t)  =  he  q  =  Z  C.e  (C7) 

"  i  =  l 

or  in  transfer  function  form,  we  have 


H(s) 


hT(sI-F)_1q  = 


N 
Z 
i  =  l 


r 

i 


(s+sj)(c>+s^*) 


;8t 


In  the  next  section,  we  discuss  three  parameter  estimation  algorithms 
applied  to  this  problem:  (1)  Prony's  technique  which  utilizes  the  models 
of  (C6),  or  (C8);  (2)  extended  Kalman  filter  technique;  and  (3)  the 
maximum  likelihood  identifier,  both  of  which  use  the  state  space  form  of 
(C6). 

3.0  Parameter  Estimation  Algorithms 

In  this  section  we  discuss  the  three  oarameter  estimation  algorithms 
employed  to  extract  the  set  of  object  parameters  {{S-. 
from  noisy  measurement  data.  The  algorithms  employed  were:  (1) 
Interactive  Prony's  technique  (If,T);  (2)  extended  Kalman  filter  (z<~); 
and  (3)  maximum  likelihood  identifier  (XXLKIO).  We  will  not  discuss  the 
mathematical  details  of  these  alqnrithms,  but  rather  include  the  primary 
references  for  the  interested  reader.  After  presenting  each  algorithm, 
we  will  compare  them  and  discuss  the  various  tradeoffs. 

Interative  Prony's  technique  (IPH  is  basically  a  linear  least  squares 
estimator  for  poles  in  the  discrete  (7  transform)  domain  [Ref.  Cl].  The 


1 


alqorithm  is  depicted  (simply)  in  Fiqure  Cl.  Dependinq  on  the  siqnal- 
to-noise  ratio  (SNR)  [Ref.  C 3 ]  either  the  impulse  resoonse  (high  SNR1  o^ 
tne  autocorrel at  ion  response  (low  SNR)  is  estimated  using  fast  Fourier 
transforms  or  sample  autocorrel at  ion  estimators,  respectively.  The 
filtered  data  is  then  "windowed",  and  ooles  estimated  from  each  data 
window  by  solving  a  set  of  linear  matrix  equations  to  obtain  linear  least 
squares  estimates  of  the  coefficients  of  a  polynomial,  the  roots  of  which 
are  the  discrete  (z  domain)  poles.  These  poles  are  then  transformed  to 
the  continuous  domain  and  identified  directly  with  the  object  resoonse 
[Ref.  Cl].  This  technique  is  repeated  by  the  orocessor  many  times  and 
"pole  clusters"  are  obtained.  It  should  be  noted  that  the  discrete  or 
samoled  data  domain  representation  is  necessary  because  of  the  use  of 
"sampled"  response  data.  Not  accountina  for  the  sampling  phenomenon, 
will  result  in  erroneous  estimates  for  the  continuous  noles. 

The  extended  Kalman  filter  (EKF)  is  basically  a  nonlinear  state 
estimation  alqorithm  which  can  be  used  to  estimate  unknown  parameters  bv 
redefining  them  as  states.  Recall  that  a  state  estimator  is  a  comouter 
algorithm  whicn  may  incorporate;  (1)  knowledqe  of  the  physical  process 
phenomonology;  (2)  knowledqe  of  the  measurement  system;  (3)  knowledge 
process  and  measurement  uncertainties  in  the  form  of  mathematical  models 
to  produce  an  estimate  of  the  state. 

'■lost  state  estimators  can  no  placed  in  a  recursive  form  with  the  various 
subtleties  emerging  in  •  o<>  '•ai cm  i -at.  mu  t  v»  current  e«,ti  y  ,  ) 
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calculations  the  new  or  corrected  state  estimate  is  ontain-'d.  :m- 
algorithm  continues  in  this  loop,  processing  measurement.  data  as  i 
oecoines  available.  Tnis  processing  is  considered  on-1  in*-  because 
be  accomplished  in  conjunction  with  t’l1  response  meas  ireaients , 
State  estimates  are  updated  in  real  time,  maon  time  i  new  a su 
becomes  available. 

The  final  a  Igor  it  if.  s  tn>  na*.  imi,-  1  1  1  i  h  and  identifier  ('-‘Oil  i. 

MXL\l!)  algorithm  is  a  co">nle<  if  f-  1  in-»  technique  whi~n  utilizes  a 
parameter  optimization  algorithm  oonft.<  around  the  rv-  to  >‘it..-iin 
parameter  estimates.  The  algorithm  ma<  imizos  ttm  1  i  •  •••!  ■  nond  r  n't 
equivalently  minimizes  the  negative  lo  - 1 ikel inonn  function  ,1'H'; 
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Figure  C2 .  ^tended  Kalman  Filter  Algorithm. 
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rf  are  the  ngr  -"“ter  ust  i1"  it*- 
0  is  a  step  parameter 
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vein  .  turn  used. 


’n'  nnpl  i  f  led  algorithm  operation  is  looii.ton  hi  :  i  :ur-.-  '  .  fh»?  -aV-  m 
i  1  t»r  is  used  to  oroducc*  -.mc-irrol  a  ted  innovations,  e ,  from  tn<? 
or  rr-1.it  on  measurements,  z.  The  1  i<ei  ihoo-.l  formation  ooeratior  ;''!!•  o 
-lieu  1 , ued  using  results  f roin  t h- •  < a  1  nan  filter.  In  some  ooti’mati  n 
1 1  aor  i  tnms  tnr  filter  is  also  used  to  calculate  eler'-'nts  in  the  .-,-einot  inn 
ritr ix,  h  (e.g.,  see  ’csl), 
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TACi.E  Cl.  IDENTIFICATION  ALGORITHM  COMPARISONS 


PT 

EKF 

MXLKID 

Problem  scope 

Linear, 

time  invariant 

Linear, 
nonlinear, 
time  varying 
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are  application,  complexity,  and  accuracy.  The  MXLKID  algorithm  is  more 
accurate;  however,  tne  price  paid  is  complexity,  computer  time  and  the 
necessity  of  running  the  algorithm  off-line,  '.e.,  with  a  complete  sjt  of 
measurements  available  beforehand. 

This  completes  our  discussion  of  the  parameter  estimation  algorithms.  In 
the  next  section  we  discuss  the  application  of  these  techniques  to  tne  :/•' 
parameteri zation  problem. 

4.0  Application  to  EM  Response  Parameter i zat ion 


The  three  analysis  methods  described  in  the  previous  sections  have  been 
applied  to  some  of  the  transient  data  described  in  Appendices  C,  '• ,  and 
E.  The  extended  Kalman  filter  and  the  maximum  likelihood  estimator  have 
been  applied  to  the  cylinder  data  only.  These  are  described  in  this 
section.  The  Prony  results  which  are  described  in  Appendix  6  include 
pole  predictions  for  all  three  objects  tested  as  part  of  this  effort. 

goth  the  extended  Kalman  filter  and  maximum  likelihood  identifier  require 
preprocess ing  of  the  data  (described  later  in  this  section)  to  obtain 
initial  estimates  for  noise  statistics,  etc.  After  comoletion  of 
preprocessing,  the  following  parameters  from  the  signal  model  (Eq.  C4) 
were  estimated;  f<r.,  w..}t  i  =  1,2,3. 

Three  sets  of  cylinder  data  were  used  in  tnis  analysis:  1)  cylinder  in 
free  space;  2)  cylinder  50  cm  from  a  qround  plane;  and  3)  cylinder  10  cm 
from  a  qround  plane.  Figure  C4  shows  the  data  records  used  which  consist 
of  the  first  20  ns  of  the  response  for  each  of  the  above  three  cases. 

•  O 

All  three  correspond  to  the  B  probe  resoonse  at  0  =  ISO  . 

The  parametric  model  for  this  experiment  with  measurement  noise  is  the 
fol lowing; 


N  s .  t 

z(  t)  =  Z  Te[C  .e  1  ]  +  v ( t ) 
i=  1  1 


ri3) 
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whe>  e 


z(t)  is  the  measured  siqnal  if  interest 
v(t)  is  the  measurement  noise 
N  =  6* 

C-  is  the  complex  residue; 

C  -  <r  a  :  +  j  1  mr  A  •  J 
s.  is  tne  complex  frequency; 

S  -  <r.  +  i,  - 

i  i  Ju  i 

[he  ro'e  extraction  problem  then  becomes; 

[ dent i f v  ( C  . ,  w  . ,  a  .  i 
for  i  -  ] ,  z ,  •  ai  yen  {/(t)J 

An  alternative  and  more  complete  noise  model  consider-;  process  as  well  as 
measurement  noise.  Process  noise  is  correlated  (as  opposed  to  independent  or 
"white"  noise  character i st i c  of  a  measurement  probe)  and  typically  has  power 
spectral  components  within  the  siqnal  bandwidth  of  interest.  Process  noise 
sources  can  include  unmodeled  IMP  reflections,  unknown  environmental 
electromaqnetics,  and  signal  mismodel ing.  Process  noises  are  modeled  as 
driving  noises,  i.e.,  random  excitations,  and  hence  are  convolved  with  the 
impulse  response  of  the  object. 

The  following  parametric  model  includes  both  process  and  measurement  noises; 

N  s-t  t  s,( t-T) 

z  ( t )  =  £  C  .  e  +  [  w  '  t  )  e  +  v '  t )  i  Z 1 4 ) 

i  - 1  ‘  T  --()  1 


*  Wp  chose  the  first  three  harmonics  'J  =  “'> ,  m.,  a  complex  conjugate  nair  of 
poles  for  each  harmonic)  to  form  our  s  i-jn-il  undel;  however,  because  of  probe 
location,  very  little  of  t fa e  second  harmnni"  was  observahl.*,  therefore  the 
forthcoming  pole  extraction  results  will  apply  to  the  fundamental  and  third 
harmonic  only. 


b'l 


Where 


wi(t)  is  process  driving  noise 
The  pole  extraction  problem  is  now: 

Identify  {C i ,w. , a. ,Cov(  Wi )} 
for  i  =  1,2,3  given  (z(t)} 

We  now  have  the  additional  task  of  identifying  the  covariance  (noise  power)  of 
the  drivinq  noise.  To  extract  ooles  from  the  three  response  data  sets  shown 
in  Piqure  C4,  tne  Prony  technique,  extended  Kalman  filter  ( )  and  two 
versions  of  the  maximum  likelihood  identifier  were  applied.  Tne  Prony 
technique  accounts  only  for  a  measurement  noise  model,  whereas  the  -  <  F  uses 
both  a  process  ane  a  measurement  noise  model.  The  maximum  likelihood 
identifier  was  designed  for  both  cases.  Version  A  accounted  only  for 
measurement  noise  whereas  Version  3  included  both  process  and  measurement 
noise  models. 

The  raw  incident  pulse  response  data  was  preprocessed  in  two  steps:  a)  low 
pass  filtering,  to  eliminate  noise  components  in  the  hiqh  frequency  bands  and 
in  which  no  signal  of  interest  lies,*  b)  discard  data  which  are  obviously 
before  and  during  the  incident  pulse  transient  (-.first  5  ns).  The  low  pass 
filter  cut-off  was  chosen  to  reflect  the  need  for  a  samplinq  rate  of  at  least 
twice  the  highest  frequency  component  of  interest  (which  is  three  times  the 
fundamental  frequency  of  -140  MHz,  i.e.,  -  420  MHz).  The  identification 
algorithms  (EKF,  MXLKIO)  were  applied  t.o  the  resultant  preprocessed  data. 

The  results  from  EKF  and  maximum  likelihood  techniques  are  shown  in  Table  7  in 
the  main  body  of  the  report.  For  a  discussion  of  these  results,  the  reader  is 
referred  to  Section  8  of  this  report. 


*  An  8-point  window  averaging  filter  was  used  (on  the  time  domain  data)  to 
produce  a  low  pass  effect. 
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Figure  C4.  FMP  response  data. 
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Appendix  E 

Crossed  Cylinder  Data 

The  measurements  were  performed  for  both  a  free-space  configuration  and  in  trie 
presence  of  a  perfectly  conducting  ground  plane.  For  the  crossed  cylinder, 
only  a  0  ano  180°  incident  field  was  used.  Figures  El  through  E4  show  the  E 
ana  6  sensor  outputs  for  the  free-space  measurements  at  0°  and  l£u°  incidence 
angles,  respectively. 

The  f requency -domain  transfer  functions  were  obtained  in  a  manner  similar  to 
that  used  for  the  cylinder  measurements;  i.e.,  the  Fourier  transforms  of  the  u 
and  6  responses  were  divided  by  the  Fourier  transform  of  the  incident  EM  pulse 
measured  at  the  base  of  the  model,  but  with  the  model  removed.  These  transfer 
functions  are  shown  in  Figures  E5  through  E8. 

The  measurements  for  the  crossed  cylinder  near  a  perfectly  conducting  ground 
plane  are  shown  first  in  the  time  domain  in  Figures  E9  through  E24,  and  then 
in  the  frequency  domain  in  Figures  E25  through  E40.  For  the  time  domain 
waveform,  only  the  first  20  ns  of  the  response  is  shown. 
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Figure  £13.  Time  domain  waveform  for  5  sensor  on  crossed  cylinder  10  cm  from 
ground  plane  (8  =  180°). 
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Fiqure  E16.  Time  domain  0  sensor  output  for  crossed  cylinder  2  m  from  qround 
plane  (e  -  180°). 
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Figure  E17.  Time  domain  waveform  for  8  sensor  on  crossed  cylinder  located 
10  cm  from  ground  plane  (9  =  0°). 
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Figure  E18.  Time  domain  waveform  for  $  sensor  on  crossed  cylinder  SO  cm  from 
gound  plane  (9  =  0°). 
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Figure  £21.  Tine  domain  waveforms  for  IT  sensor  on  crossed  cylinder  located 
10  cm  from  ground  plane  (9  -  130°). 


10  1  ‘  -  -  *  1  J  1  '  J  -l  -  1 

0  4  8  12  16  20 


lime  -  ns 

Figure  £22.  Time  domain  waveform  for  IT  sensor  on  crossed  cylinder  located 
50  cm  from  ground  olan^  (B  -  180°). 
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Figure  £23.  Time  contain  waveform  for  3  sensor  on  crossed  cylinder  located  1  m 
from  ground  plane  (9  =  110°). 
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Figure  £24.  Time  domain  waveform  for  3  sensor  or,  crossed  cvl  i  .  v->r  iocated  '  m 

from  ground  olane  (9  = 
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Figure  E2S.  Magnitude  of  E/E^  for  crossed  cylinder  17  c-n  fr<vn  a  perfect 

ground  plane  (6  =  0°) . 
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Figure  £27.  Magnitude  nf  * /*  .  for  crossed  cvlinder  '  m  frr  r  -|**ound  olane  1 

(0  -  0°). 
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Appendix  : 

74/  Sralr- Model  A  i  rv  r  u  t  0,5ta 

Two  types  of  measurements  were  made  mi  tn--  air.. aft  model .  The  first  was  for 
a  free-soace  conf igurat ion  with  th>  in  i r  ,*  i.-rtr  i<:  field  parallel  to  the 
wings  from  the  top.  Charge  neasuomumts  were  made  at  the  top  of  the  wina  tin 
and  midway  along  the  wing,  .voile  the  sort  <ci;  .a  rent  was  measured  at  the  same 
midwing  point  and  at  tne  winp  root.  The  siirfac-  currents  were  also  measured 
at  tne  test  point  on  tne  fuselage. 

The  other  half  =  ■  *  the  / 4 /  m.d*  i  -ms  mount  <■  d  a  ,  before.  However,  this  t  i<-<*  the 
sensor  location,  wr.*  •  ar  :(**d  on  the  pot  tom  of  tie1  wing.  The  wheel  wells  were 
again  closed  anr  'in"i-  r,  mi  the  monel  was  placed  an  equivalent  distance 
from  tne  ground  >lan>  is  r  it  were  sitting  on  its  wheels.  This  olaced  tne 
bottom  of  the  fiselagc  a*  rim  wing-.  ?  cm  from  the  ground  plane.  With  this 
conf iqurat ion ,  it  was  necessary  to  run  th>-  coaxial  cable  from  the  0  sensor  at 
the  wing  tin  along  the  ton  (incident  field  side)  of  the  wing.  We  used  a  very 
thin  section  of  tne  semiri  iid  coaxial  cable  to  connect  the  sensor,  and  covered 
the  coaxial  with  the  conductive  cooper  tap.-.  Wo  noticeable  change  was 
observed  if  the  coaxial  section  was  covered  or  not. 

The  time-domain  results  of  the  free-soace  measurements  are  shown  first  iri 
Figures  FI  through  F5  for  each  of  the  test  points.  As  before,  the  temporal 
measurements  were  then  transformed  to  the  frequency  domain  to  obtain  tne 
desired  frequency-domain  functions.  These  results  are  shown  in  Figures  FA 

through  F10. 

Figures  FIT  and  -XL  show  the  time-domain  responses  obtained  for  two  specified 
test  points  with  the  model  located  over  the  perfect  ground.  ror  comparison, 
we  also  measured  the  mr.de 1  response  with  the  ground  plane  removed.  These 
waveforms  are  shown  in  cigures  F13  and  FIT.  (>n  all  cases,  only  the  first  90 
ns  of  the  response  are  shown.) 


Finally,  the  time-domain  waveforms  were  converted  to  the  frequency  domain  and 
normalized  to  obtain  the  transfer  functions,  shown  in  Fjqures  c15  through  -13, 
which  correspond  to  the  ground  plane  and  uottom  side  free-space  results. 


figure  F3.  Transient.  response  rtf  •>  senso*-  at  t -in  nidwing  for  model 
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figure  F4.  Transient  reeno*  so  nr  sons  r  .it  -'inn  '•ait.  for  «•«•%  j 
spare. 
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Figure  F6.  Magnitude  of 
free-space  configuration. 
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Figure  Fl<>.  Transient  response  of  1  on  nottom  surf  ace  of  wine  wife  aircraft 
over  a  perfectly  conducting  ground  plan. 
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Fiqure  Fli.  Transient  response  of  0  sensor  on  not  ton  tin  of  sjnq  witn 
aircraft  in  free-space  mode. 
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-iquro  ‘•'14.  Tr  ansient  response  of  8  sensor  ir.  oott.rsi  surf  y  ,■  of  w  mn 
aircraft  in  frpe-snar.o  'node. 
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PRONY  PROCESS  I  Nil  Oc  rut  RAN3[rNT  OArA 

Tne  Singularity  Expansion  Method,  (Ref.  6 !  and  ■  ?),  and  the  us.?  of  Pronv's 
algorithm  (Ref.  if)  to  extract  -sinqulartit  i<>s  tram  transient  time-domain  data 
has  seen  considerable  attention  in  electromagnetics  durinq  the  past  3-4  years 
(Refs.  G4  and  G5).  Basically,  Prony's  alqontivn  is  a  method  of  determinim 
tne  coefficients  a-  and  S-  in  the  c omn I e *  exponential  series  representation 
of  a  time  series: 


f(t  .1  -  ^  Ro?sit.,  j  I, - m  1  time  sample  (51  I 

i  =  l 

where  the  S(  =  tr  +  ju  are  the  complex  singularities  or  noles  of  the  time 
function,  and  the  R.  are  the  complex  residues.  The  s  - ' s  otitained  from 
Prony's  procedure  have  particular  significance  in  electromagnetics  because 
they  are  characteristic  of  the  natural  frequencies  of  oscillation  found  for 
bodies,  such  as  the  cylinder  and  the  aircraft,  wnen  they  are  excited  by  an 
impulse  of  energy. 

Once  the  complex  £r.,  s.J  pole-residue  set  is  obtained,  the  corresponding 
frequency-domain  response  is  given  by  the  Laplace  transform  of  the  above  time 
series  representation: 

F(S1  ^  j -  V  •  (G’l 

Several  computer  programs  have  been  writt'-n  t.o  dn  the  Prony  Processing  of 
time-domain  data,  including  the  SrMP‘ X  code  used  at  I.LNL  (Ref.  66) .  All  known 
3rony  techniques,  however,  suffer  from  deuraded  performance  if  noise  is 
present  in  the  input  da:a,  part icularl y  if  tfn>  data  is  indiscriminately 
processed.  At  LLNL,  work  has  been  done  on  lindinq  ways  to  minimize  the 
detrimental  effects  of  the  noise  on  the  Prony  process.  These  techniques 
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include  the  use  of  prefilters  to  limit  unwonted  high-frequency  noise  content 
in  tne  waveforms,  and  the  use  of  multiple-run  averaging  to  smooth  out  the 
effects  of  noise. 

The  latest  program  development  has  been  incorporation  of  the  Prony  orocessor 
into  a  large  general  purpose  data  analysis  program  called  GPDAP  (Ref.  G7). 
GPDA3  runs  on  a  committed  mini  computer  system  and  has  utilities  such  as  fast 
Fourier  transforms  and  filtering  urograms  which  allow  the  user  to  easily  and 
quickly  manipulate  the  process  data  and  then  plot  the  results  on-line.  The 
flexibility  of  tnis  proqrant,  coupled  witn  a  high  degree  of  interactivity, 
gives  the  user  new  insiqht  into  the  application  of  the  Prony  program  to 
process  real  experimental  data,  because  the  user  can  quickly  see  the  outcome 
of  his  work  and  make  any  desired  changes  while  the  many  oaramenters  are  fresh 
in  his  mind.  The  results  presented  in  the  next  subsections  are  those  obtained 
with  this  new  system. 
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Prony  Analysis  for  the  Cylinder  Model 


Experimental  measurements  on  the  cylinder,  crossed  cylinder,  and  the 
scale-model  747  aircraft  were  all  analyzed  with  Prony's  method.  In  this 
subsection  the  results  for  the  cylinder  in  free  space  and  near  a  perfect 
qround  are  oresented. 

All  of  the  measurements  we  performed  were  made  on  a  plane  of  symmetry,  wbicn 
implies  that  only  odd  modes  are  excited  by  the  incident  field,  as  illustrated 
in  Figure  Gl.  I  should  he  noted  from  this  fiqure  that  some  modes  may  not  have 
large  modal  amplitudes  at  the  point  of  observation  used  in  these  experiments; 
and,  in  fact,  some  of  the  modes  may  be  unobservable  or  at  least  below  the 
useful  si gnal -to-noise  ratio  needed  to  obtain  qood  results  from  the  Prony 

processor.  The  transient  response  waveforms  used  for  these  brony  calculations 

•  • 

were  the  direct  3  and  0  sensor  outputs,  and  as  such  they  were  not  normalized 
to  the  incident  exciting  field. 

While  all  of  the  measurements  obtained  from  this  series  of  experiments  could 
be  analyzed  by  the  Prony  program,  we  have  chosen  a  subset  of  the  measurements 
for  the  response  near  a  perfect  ground  plane  to  determine  the  effects  of  tne 

nearby  reflector.  Additionally,  measurements  for  the  cases  of  180°  incidence 

•  • 

were  selected  for  the  3  and  0  sensors.  ciqures  G2  (al,  (b),  and  (cl 

■ 

illustrate  the  transient  waveforms  from  the  3  sensor  with  the  cylinder  located 
10  cm  ( 2 a  1 ,  50  cm  (10a)  and  2  m  (40a)  away  from  the  ground  plane,  respectively 

The  choice  of  the  time-window  to  use  for  the  P  ony  analysis  depends  upon 
several  factors.  First,  because,  the  primary  interest  is  finding  the  natural 
modes  of  the  body,  we  select  a  f ree-resoonse  or  undriven  portion  of  the 
waveform  after  the  incident  pulse  has  passed,  second,  to  observe  the  effect 
of  the  ground  plane,  we  must  wait  a  time  sufficient  for  the  incident  pulse  to 
strike  the  body,  reach  the  ground  plane,  and  then  reflect  back  to  the  body 
again.  (This  ground  plane  reflection  can  clearly  be  seen  occurring  at 
different  times  in  Figure  G2.)  Finally,  the  duration  of  the  measurment 
interval  must  be  sufficiently  lonq  to  allow  a  separation  of  r.losplv  spaced 
modes,  the  window  of  observation  must  be  ?1/Af. 


These  requirements  mean  that  we  can  determine  the  effects  of  the  ground  plane 
for  the  cases  of  the  object  at  10  cm  and  50  cm  spacing  of  the  qround  plane, 
but  are  unable  to  analyze  the  cases  of  1  m  and  ?  m  because  the  measurement 
interval  was  insufficiently  long  to  meet  the  observation  period  criteria. 
Instead,  the  case  of  the  2  m  spacing  is  used  as  a  free-space  case  for  the 
period  between  the  initial  incident  pulse  and  the  reflected  pulse  from  the 
ground  plane. 

To  condition  the  time-domain  waveform  for  nrony  processing,  the  Fourier 
transform  of  the  signal  if  first  obtained.  The  spectrum  is  then  truncated 
above  a  chosen  cutoff  frequency  ( ? - 08  GHz).  Next,  an  inverse  rourier 
transform  is  made  of  the  truncated  data  to  get  the  siqnal  back  to  a  filtered 
time-domain  waveform  which  is  then  used  as  the  Prony  input. 

The  Prony  processor  is  then  used  on  the  filtered  data  to  form  a  complex 
exponential  curve  fit  to  the  time  series  data.  This  procedure  determines 
first  o.'s,  which  are  the  complex  poles  of  the  structure  in  conjugate 
pairs.  With  these  poles  and  the  original  data,  the  complex  residues,  S., 
are  computed  to  fit  the  time  data.  The  primary  interest  here  is  in  the 
location  of  the  poles  in  the  complex  plane  { s  =  <j  +  jw)  .  cor  the  cylinder,  we 
show  in  Figures  G3  (a)  through  (c)  the  location  of  the  resulting  poles  in  the 
upper  left-half  plane  determined  for  the  cases  of  the  qround  plane  10  cm,  50 
cm  away  from  the  cylinder,  and  also  the  free  space  result.  The  pole 
locations  result  from  a  curve-f i tt i nq  process  where  we  intially  asked  for  30 
pole  terms  to  fit  the  data,  and  then  selected  closely  grouped  clusters  of 
significant  ooles  by  drawinq  boxes  around  the  clusters,  as  shown  in  Figure 
G4.  Here  each  dot  on  the  plot  represent  a  pole  location,  and  the  cluster 
result  from  10  separate  curve  fits  to  the  time  data  usinq  a  si idinq-window 
approach.  All  of  the  poles  that  fall  within  the  boxes  are  averaged  to  qive  a 
mean  pole  location.  These  mean  pole  locations  are  indicated  by  the  dots  in 
Figures  G3  (a)  through  (c),  and  the  boxes  drawn  around  the  pol°s  are  the 
variances  of  the  poles  in  both  frequency  and  damping.  Note  that  generally  the 
variance  in  the  value  of  the  damping  coefficient  is  much  larger  than  the 
variance  in  the  complex  frequency  of  tne  pole. 
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Once  the  significant  poles  are  selected,  they  dire  used  in  the  curvefittinq 
process  to  compute  the  residues  (S^.  The  values  of  the  poles  and  residues 
for  the  three  cases  of  the  3  sensor  on  the  cylinder  are  tabulated  in  Table  Gl. 

Figure  G5  shows  a  composite  plot  of  the  poles  for  the  three  cases  analyzed. 

The  two  lowest-order  oole  sets  are  well  behaved,  while  the  results  for 
frequencies  above  500  MHz  are  quite  scattered.  The  trend,  however,  can  be 
seen  as  the  cylinder  nears  the  ground  plane. 

Once  the  poles  and  residues  are  found,  they  may  be  used  in  the  complex 
exponential  series  representation  of  the  time  waveform.  Me  can  see  how  well 
the  curve  fit  preforms  by  examining  the  plots  shown  in  Figures  G6  (a)  through 
(c),  which  show  both  the  original  data  (solid  curve)  and  the  curve  computed 
form  the  poles  and  residues  (dotted  curve). 


•  • 

In  a  manner  very  similar  to  that  used  for  the  B  measurments,  the  D  measure¬ 
ments  on  the  cylinder  were  also  analyzed  with  Prony.  The  same  cases  were  run, 
and  the  reader  can  refer  back  to  the  measurements  section  on  the  cylinder  to 
see  the  transient  waveforms  used  here.  Again,  30  poles  were  initially 
requested  in  the  curve-fitting  process,  and  10  sample  data  were  made  to  form 
pole  clusters.  The  resulting  pole  plots  are  shown  in  Figures  G7  (a)  through 
(c)  for  the  two  ground-plane  soacings  and  the  free  space  case.  The  poles  and 
residues  are  also  tabulated  in  Table  G2  for  the  0  sensor  on  the  cylinder, 
along  with  a  composite  pole  plot  for  the  three  cases  shown  in  Figure  G8.  To 
see  now  well  the  selected  poles  fit  the  input  waveform,  the  reconstructed 
curves  are  shown  along  with  the  input  in  Fiqures  G9  (a)  through  (c). 

Prony  Analysis  for  the  Crossed  Cylinder  Model 

[n  addition  to  the  cylinder,  Prony's  method  was  also  used  on  the  experimental 
measurements  obtained  for  the  crossed  cylinder.  Aqain  a  selction  of  the  runs 
in  the  presence  of  a  perfectly  conducting  ground  plane  was  made.  The  cases  of 
0  -  180  incidence  with  the  B  and  0  sensors  were  also  selected,  as  with  the 
cylinder.  Fiqures  G10  (a)  throuqh  (c)  show  the  transient  B  waveforms,  while 
Figures  Gl 1  (a)  throuqh  (c)  show  the  analogous  measurements  for  the  0  sensor. 
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TABLE  G2 . 

POLE  VALUES  FOo  CYLINDER  D-SENSOO  (&  =  180°). 
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As  with  the  cylinder  measurements,  the  Fourier  transforms  of  the  time  data 
were  truncated  above  2  GHz  and  then  inverse  transformed  to  orovide  filtered 
time-domain  data  for  the  Prony  proqram.  The  results  of  the  prony  analysis  are 
presented  next. 

For  each  of  the  3  measurements,  a  30-pole  curve  fit  was  selected,  and  the 
results  for  a  10-run  sliding  window  were  plotted  on  one  pole  nlot  to  indicate 
clustering  if  the  pole  locations.  The  cases  where  the  ground  olane  was 
located  10  cm  and  50  cm  away  from  the  model  were  used  as  ground-plane 
measurements,  while  the  case  where  the  ground  olane  was  located  ?  m  away 
the  model  was  used  as  the  free-space  measurement  because  relection  back  fm 
the  ground  plane  did  not  occur  during  the  data  time  interval  used  For  the 
Prony  processor.  Fiqures  Gl?  (a)  thrpuqh  (c)  show  the  oole  plots  alonq  with 
the  variance  in  the  pole  locations  for  the  three  cases.  Table  G3  is  a 
tabulation  of  the  mean  oole  values  along  with  their  residues.  These  poles  a1'-? 
all  plotted  in  Figure  G13.  The  reconstructed  P  waveforms  using  the  selected 
noles  are  shown  in  Fiqures  G 1 4  (a)  throuqh  (c). 

Analysis  for  the  n  sensor  output  on  the  crossed  cylinder  follow  the  same  oath 
as  those  before.  Figure  Gin  (a)  throuqh  (cl  show  the  pole  plots  for  the 
individual  measurements.  Table  G4  summarizes  the  mean  pole  values  found, 
while  the  composite  plot  of  Figure  GI6  shows  the  poles  for  all  three  runs. 

The  reconstructed  waveforms  using  the  poles  and  residues  are  shown  in  Figures 
G17  (a)  through  (c). 

Prony  Analysis  of  the  747  Aircraf*  Model  Measurements 

Included  in  the  time-domain  measurements  were  two  exneriments  nerformed  on  a 
1:100  scale  model  of  a  747  aircraft.  The  ronf iuurat ion  of  the  experiments  is 
shown  in  Fiqure  G1S  where  we  measured  the  surface  charge  density  on  the  bottom 
of  the  winq  tip  and  the  axial  surface  current  density  on  the  bottom  of  the 
winq  midway  between  the  wing  tip  and  the  fuselage.  These  two  measurements 
were  repeated  for  the  case  where  the  aircraft  would  be  located  the  same 
distance  above  tne  qround  plane  as  it  would  witn  the  wheels  down,  and  the  case 
where-  the  ground  plane  was  removed  for  a  fme-space  measurement  .  In  both 
cases  the  wheels  were  up  and  thp  wheelwell  doors  were  closed  and  painted  over 
with  silver  conductive  oaint. 
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TABLE  03. 

POLES  AND  RESI'JDES  FOB  CROSSED  CYLINDER  3  SENSOR  (9 
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TABLE  G4. 

CROSSED  CYLINDER  POLES  AND  RESIOiJES  0  SENSOP  (9 


The  Prony  poles  for  the  surface  charge  density  on  the  wing  tip  were  determine:) 
first.  Forty  data  points  between  0.3b  and  1^.7  ps  were  used  as  input  to 
Prony,  and  the  process  was  repeated  over  1 -j  data  windows  as  with  the  cylinder 
and  crossed  cylinder.  The  resulting  mu i t ip le-pole  plots  for  the  upper 
left-half  plane  are  shown  in  ^iqums  ..I'd  .'<•!  and  (h)  for  the  cases  of  the 
aircraft  on  the  ground  and  in  free  space.  rh.  boxes  are  shown  drawn  around 
the  most  dominant  pole  clusters,  and  the  m.-an  values  of  the  boxed  poles  are 
used  to  calculate  the  residues.  These  results  are  tabulated  in  Table  GS.  The 
resulting  fit  to  the  sampled  data  is  shown  in  reconstructed  curves  in  Figures 
G70  (a)  and  (b),  where  the  solid  line  represents  trie  original  data  and  tie- 
dotted  curve  is  the  resulting  curve  fit. 

Finally,  the  measured  results  for  trie  d  sen  >or  on  the  mil-losition  of  the 
aircraft  wing  were  analyzed. 

The  Prony  cluster  plots  are  shown  in  F  iqur<  .i.M  ■,  a)  and  (  o)  wnere  we  have 
selected  the  two  low-order  dominant  poles.  Trie  nean  values  of  the  poles  eiil 
their  residues  are  tabulated  in  Table  Go.  To  show  how  well  these  two  poles 
the  measured  data,  the  waveform  reconstruc t ed  is  shown  last  in  Figures  G?P  (a' 
and  (b). 
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P90NY  POLES  FULL-SCALE  747  MODEL  ON  THE  GROUND  AND  IN  S»ACE 

THESE  POLES  ARE  FOR  THE  CHARGE  DENSITY  ON  THE  ROTTOM  SUReACE  WING  TIP 
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"f ree-soace"  1  configurations.! 
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riqure  Go.  (b)  Measured  and  reconstruct.:*')  wav  >'oms  for  cv  i^-;er  SO  cnr  fro- 
qroiind  plane. 
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-igure  VI.  {:)  Measured  am;  reconstruct.  •!  wav-forms  for  r  v  r  >  •>  t>f  • 
ground  plane.  (Note:  Only  free-soace  portion  of  waveform  is  used  for  tm* 
l’rony  anal  /sis. ) 
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Figure  G13.  Pole  locations  for  crossed  cylinder  3  sensor  (0 
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Figure  G 1 4 .  (a)  Measured  and  reconstructed  8  waveforns  for  crossed  cylinder 

10  cm  from  ground  plane. 
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